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Abstract 

In  previous  work  we  introduced  a method  of  using  polynomial  splines  with  appropriate 
discontinuities  to  approximate  a piecewise  smooth  function  / with  jump  discontinuities  of 
/ and  }' . The  information  used  is  location  of  discontinuities,  and  low  order,  possibly  noisy 
Fourier  coefficients.  The  number  of  discontinuities  was  limited  to  two  at  most,  and  the 
discontinuities  needed  to  lie  at  meshpoints  in  a uniform  mesh.  We  showed  that  the  linear 
operator  corresponding  to  the  method  is  1,2-bounded  with  a modest  bound,  and  thus 
that  the  method  is  Z/2-robust  in  the  presence  of  noise.  In  the  present  paper  we  develop 
a new  method  of  analysis  which  enables  us  to  determine  operator  bounds  that  are  valid 
for  arbitrarily  many  discontinuities.  The  new  analysis  allows  discontinuities  to  be  placed 
arbitrarily.  Given  a placement,  an  initially  uniform  spline  mesh  of  width  h must  be  used 
such  that  nearest  meshpoints  to  discontinuities  are  at  least  4 h apart  (discontinuities  then 
replace  these  meshpoints);  the  number  of  available  Fourier  coefficients  must  be  at  least 
three  times  the  number  of  mesh  intervals  in  a period.  The  previous  work  was  restricted 
to  quadratic  splines;  the  present  work  includes  cubic  splines.  Much  of  the  analysis  uses 
exact  computations  with  a computer  algebra  system.  We  give  an  example  to  illustrate 
the  accuracy  of  the  method  using  noisy  Fourier  coefficients. 

1 Introduction 

We  consider  approximating  a function  / when  the  information  consists  of  low  order,  pos- 
sibly noisy  Fourier  coefficients,  and  knowledge  that  / is  smooth  except  for  jumps  of  / or 
f at  known  locations  but  unknown  magnitudes.  We  will  work  with  a method,  introduced 
in  [10],  which  amounts  to  linear  least  squares  fitting  of  the  available  coefficients  with  the 
coefficients  of  splines  with  appropriately  placed  discontinuities.  Since  we  anticipate  ap- 
plications to  ill- posed  problems  where  boundedness  of  the  solution  operator  is  crucial, 
we  develop  a method  for  bounding  the  norm  of  this  operator.  The  bounding  method 
depends  heavily  on  exact  computations  in  certain  spline  spaces.  These  computations  are 
fundamentally  finite  dimensional  linear  algebra  with  rational  integer  coefficients.  Their 
goal  is  to  develop  upper  bounds  for  the  norms  of  certain  projector  operators  whose 
norms  are  naturally  expressed  in  terms  of  generalized  eigenvalues,  and  to  prove  by  exact 
computation  that  the  bounds  are  correct.  A computer  algebra  system  is  used  for  the 
computations.  The  programming  is  detailed  in  [9]. 
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In  [10]  we  obtained  bounds  under  much  more  restrictive  conditions  than  in  the  present 
paper.  In  [10]  the  splines  were  quadratic  only,  while  here  results  also  are  given  for  cu- 
bic splines.  The  analysis  in  [10]  required  all  knots  of  the  approximating  splines  to  be 
uniformly  spaced,  and  since  the  discontinuities  are  at  the  knots,  the  location  of  discon- 
tinuities was  limited.  Further,  in  [10]  the  estimation  process  is  linear  in  the  total  number 
of  discontinuities,  and  produces  results  unacceptably  large  for  cases  with  more  than  one 
discontinuity  of  / and  two  of  /'. 

Others  ([2,  3,  4,  5])  have  addressed  questions  of  accurate  approximations  to  functions 
with  discontinuities  given  Fourier  coefficients  as  information.  In  [8]  we  give  examples 
which  show  that  those  methods  can  substantially  magnify  noise  in  the  coefficients;  our 
main  concern  here  is  to  prove  robustness  of  our  method.  We  illustrate  with  an  example 
in  Section  5. 

2 General  linear  space-theoretic  results 

Let  V be  a real  Hilbert  space  with  inner  product  { , ).  We  will  denote  the  norm 

associated  with  ( , ) by  ||  ||.  Let  V and  Q be  closed  subspaces  of  V;  suppose  P is  the 
orthogonal  projector  on  V.  Here,  as  in  [10],  we  deal  with  the  approximation  /*  obtained 
as  the  solution  to  the  constrained  least  squares  problem 

min  ||P/*  -P/H,/*  G Q. 

Assuming  that  P is  invertible  as  a mapping  on  Q,  we  denote  by  P+  the  mapping  from 
P(Q)  to  Q which  inverts  P.  It  is  not  hard  to  verify  that  /*  = P+RPf  where  R is  the 
orthogonal  projector  on  P(Q).  Let  A denote  the  operator  that  takes  / to  /*. 

Theorem  2.1  Let  C be  a mapping  from  V to  Q.  Let  e be  T-periodic  and  in  L2(0,  T). 
Then 


II A(Pf  + e)  - /||  < (IIP+H  + 1)|| Cf  - /||  + IIP+H  ||e||. 

Proof:  A(Pf  + e)  = Af  + Ae.  \\Af  - /||  < ||A/  - C/||  + ||Ctf  - /||  = || A(f  -Cf) || 
+ ||C-/  - /||  < (||A||  + 1)||/  - C/||.  ||A||  = ||P+PP||  < IIP+H  because  P and  R are 
orthogonal  projections.  □ 

A main  objective  of  the  following  work  will  be  to  bound  ||P+||.  This  will  be  done 
by  establishing  upper  bounds  for  ||7  — P||  as  a mapping  on  Q.  From  these,  bounds  can 
easily  be  derived  for  ||P+||. 

Theorem  2.2  Let  77  < 1 exist  such  that  ||(7  - P)qj|  < r;||g||,  for  all  q € Q.  Then  P is 
injective  as  a mapping  on  Q and  for  all  h € P(Q),  P+,  the  inverse  of  the  restriction  of 
P to  Q,  satisfies 

IIP+hll^^Hhll2. 
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We  will  obtain  bounds  for  ||/  - P||  by  considering  the  projector  perpendicular  to  a 
spline  space  Q which  is  more  tractable  than  PV,  and  on  which  7 - P is  small.  In  the  next 
section,  Q is  the  approximating  spline  space,  S a subspace  of  maximally  continuous 
splines,  and  Q is  a space  of  maximally  continuous  splines  whose  knots  are  in  a mesh 
refining  the  mesh  for  the  members  of  <S.  S and  Q have  orthogonal  projectors  S and  G, 
respectively.  The  following  estimates  ||/  — P||  in  terms  of  || I - G||. 

Theorem  2.3  Suppose  ||(/-P)g||  < 7?o!i<?||  for  all g £ Q.  Suppose  ||(/-G)q,||  < r?i||<?|| 
for  all  q £ Q.  Then  ||(7  - P)g| | < (g0  + WIMI  for  for  <?  £ 2- 

Proof:  For  q £ Q,  ||(J  - P)q\\  < ||(7  - P)G?||  + ||(7  - P)(7  - G)q\\. 

||(7  - P)Gq\\  < Vo\\Gq\\  < r?0|M|,  and  ||(7  - P)(7  - G)q\\  < ||(7  - G)q\\  < m\\q\\.  □ 


Theorem  2.4  enables  us  to  bound  ||7— G||  on  Q by  instead  bounding  projectors  orthogonal 
to  small  subspaces  of  Q,  restricted  to  small  subspaces  of  Q. 


Theorem  2.4  Let  Q and  S be  closed  subspaces  ofV  with  S C Q n Q.  Let  Vi,  V2, . . . , Vr 
be  nonzero  mutually  orthogonal  subspaces  ofV.  Let  Qi  Q QC\Vi,  1 < i < r be  nonzero 

closed  subspaces  such  that  Q = S + Qi  + 0.2-\ 1 - Qr-  Let  Qi  C Q n Hi  C fl  Vj, 

1 < i < r be  nonzero  closed  subspaces  with  orthogonal  projectors  Gi,Hi.  Let  v be  a con- 
stant such  that  ||(7  - Gi)qj]\ 2 < i/||77j<7,:||2  for  all  qi  £ Qi,  1 <i<r.  Then  ||(7  - G)q\\2  < 
u\\q\\2  for  all  q € Q. 


Proof:  q £ Q can  be  written  q = s + v where  s £ S and  v = <71  + 92  H \-qT,  qi  £ Qi, 

1 < i < r.  ||(7  - G)q\\  — ||(7  - G)n||  since  S C Q.  Let  F = G\  + G2  + ■ • • + Gr.  Since 
Gx  + g2  + ■ ■ • + Qr  c g,  ||(7  - GHI2  < 11(7  - P)n||2  = E[=1  ll(/  - Gi)q,\\2,  the  latter 
equality  because  of  orthogonality  of  the  Qi.  ||g||2  > ||(7  - S)v\\2  > ||  ]C[=i  HM\2  — 
E[=il|77igi||2,  since  ^ an<i  the  Hi  are  orthogonal.  If  all  77^  = 0 the 

hypothesis  implies  all  (7  - Gi)qi  = 0.  The  above  then  implies  (7  - G)q  = 0,  and  the 
conclusion  is  true.  We  proceed  assuming  H,q,  ^ 0 for  some  i and  let  J\f  be  the  set  of  all 
those  i.  Then 

||(J-G)<?||2  < 1K7~  Gj)gi\\ 


An  elementary  argument  shows  the  quotient  of  sums  is  < v since  for  each  i £ Af, 

11(7  — Gi)qi\\2 /\\Hiqi\\2  < v.  □ 


3 Bounds  for  restricted  projectors 

Below,  we  specialize  the  spaces  of  the  last  section,  and  get  our  main  results.  Let  T > 0 
be  a fixed  period.  We  take  V to  be  the  space  of  real-valued  T-periodic  functions  which 
belong  7,2(7)  for  some,  and  thus  every,  period  interval  7.  On  V and  its  subspaces  we  define 
the  inner  product  ( f,g ) = fjf(t)g(t)dt,  I a period  interval.  The  other  realizations  are 
defined  in  the  statements  and  proofs  of  the  following  results.  Lemma  3.1  sets  up  an 
application  of  Theorem  2.4;  Theorem  3.2  uses  this,  together  with  Theorem  2.2,  to  get 
our  main  result. 

Lemma  3.1  Let  X be  a finite  set  of  points  in  [0 ,T).  Let  N > 4 be  an  integer.  Let 
K = {iT/N,  0 < i < N}:  for  each  x £ X,  let  kx  be  a member  of  K closest  to  x where 
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0 is  identified  with  T.  Assume  N large  enough  that  between  any  two  distinct  kx  are  at 
least  three  other  members  of  K.  Let  Kx  result  from  substituting  in  K each  x £ X for  its 
kx.  For  m = 3,4  let  Q be  the  space  ofm-th  order  T -periodic  polynomial  splines  with  Kx 
as  knots  and  with  continuity  Cm~2  at  all  knots  except  the  x £ X,  where  no  continuity 
is  required.  Let  Q be  the  space  of  m-th  order  periodic  splines  with  knots  in  [0,  T ) at  the 
points  {IT/ (37V),  0 < i < 37V},  and  let  G be  the  orthogonal  projector  on  Q.  Then  I — G 
restricted  to  Q satisfies  ||/  — G\\2  < .69  ifm  = 3,  and  ||/  — GUI  < .9  if  m — 4. 

Proof:  Let  S be  the  subspace  of  Q consisting  of  those  splines  which  are  C°°  at  the 
kx.  Clearly  S C Q.  Let  h = T/N.  Fix  x = X{  £ X = {x\,  x2,  ■ ■ . , xr } and  let  y0  = Xi, 
ya  — kXi  — ah,  a = —2, —1,1,2.  Take  Vl  to  be  the  subspace  of  V consisting  of  those 
functions  with  support  in  [y-2,y2]  and  its  T- translates. 

For  m = 3 let  ji  and  j2  be  B-splines  with  knots  y~i,yo,Vo,yo  and  yo,yo,yo,yT,  let 
js  be  the  difference  of  the  B-splines  with  knots  y-2,y-i,yo,  2/i  and  2/-i,J/o>J/i>2/2  (see  [1] 
for  explanation  of  multiplicity  versus  degree  of  continuity).  For  m = 4 let  j\  and  j2  be 
B-splines  with  knots  y~i,yo,yo,Vo,yo  and  yo,yo,yo,Vo,yT>  let  h be  the  difference  of  the 
B-splines  with  knots  2/-2)2/-i,2/o,2/o,?/i  and  y-i,yo,yo,yi,y2'>  and  let  j'4  be  the  B-spline  with 
knots  2/_2,j/- 1,2/0)  y 1,2/2-  Since  y2~y~2  <T  we  may  identify  the  ja  with  their  T-periodic 
extensions. 

Let  Qi  be  the  space  of  splines  whose  generic  member  is  qi  = )T"'=1  caja  for  constants 
ca.  For  each  i,  nonzero  members  of  Qt  have  continuity  from  Cm~2  through  full  discon- 
tinuity at  Xi,  while  members  of  <S  are  C°°  at  a It  follows  that  <Sn(Qi  + Q2-l hQr)  = 0 

and  Q.  = S + Qi  + • • ■ + Qr- 

Let  Qi  be  the  subspace  of  G with  basis  the  Cm~2  periodic  B-splines  whose  knots 
in  the  period  containing  [y^2,y2]  are  length  m+  1 sublists  of  consecutive  knots  from 
the  list  ( ah/3  + kx,  —6  < a < 6).  Let  Hi  be  the  space  of  those  m-th  order  periodic 
splines  which  in  [— T/2  + kx,T/2  + kx\  have  support  in  [y^2,y2],  which  have  knots  at 
the  yi,  i ^ 0 and  at  x,  are  Cm~2  at  y_  1 and  y\,  which  may  be  fully  discontinuous  at 
y-2,V2,  and  x,  and  which  are  orthogonal  to  all  members  of  S.  \\(I  — Gi)qi\\2 /\\Hiqi\\2  is  a 
ratio  of  quadratic  forms  in  the  ca.  An  upper  bound  v for  it  can  be  obtained  as  an  upper 
bound  for  the  eigenvalues  of  the  pencil  A - A B where  aap  = ((/  — Gi)ja,  (I  — Gi)jp), 
bap  = ( Hija,Hijp ),  1 < a,/3,  < m. 

In  [9]  explicit  bases  for  the  spaces  Qi  and  Hi  are  calculated  as  m-th  order  splines. 
From  their  definitions  ([1]),  B-splines  are  rational  functions  of  the  knots,  and  thus  are 
also  inner  products  of  B-splines.  The  null-basis  and  orthogonal  projection  calculations 
in  [9]  use  standard  methods  which  involve  only  rational  operations.  Thus  the  ( I — Gi)ja 
and  Hija  and  then  the  aap  and  bap  are  rational  functions  of  the  knots  of  qi,  so  long 
as  x remains  in  [ kx,kx  + h/ 3].  When  x crosses  into  [kx  + h/3,  kx  + h/2],  thus  crossing 
knots  for  splines  in  Gi,  the  rational  functions  change,  so  in  general  the  matrix  entries 
are  piecewise  rational  functions  of  x. 

Let  v be  a conjectured  upper  bound  for  the  maximum  eigenvalue  Amax  of  A — A B (in 
[9]  a floating  point  approximation  to  Xmax  is  plotted  as  a function  of  a;;  u is  determined 
from  inspecting  this  plot).  For  computational  convenience  in  [9]  we  represent  x as  2eh/3+ 
kx,  0 < e < 1/2  for  x < kxFh/ 3,  and  as  (1  + e)h/3+  kx,  0 < e < 1/2  for  kx  + h/ 3 < x < 
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kx  + /i/2.  For  further  convenience  we  take  kx  = 0,  clearly  losing  no  generality.  We  have 
represented  only  x > kx,  but  because  of  symmetry,  x < kx  produces  the  same  bounds. 

Since  h is  a linear  factor  in  all  knots  in  the  calculation,  we  see  that  aap  and  bap 
can  be  written  as  h multiplying  piecewise  rational  functions  of  e (with  integer  rational 
coefficients).  The  determinant  of  A - vB  is  thus  hm  times  a piecewise  rational  function 
of  e.  The  MAXRAT  algorithm  ([9])  proves  that  its  reciprocal  is  bounded  as  a function 
of  e in  the  appropriate  ranges,  so  the  determinant  itself  is  bounded  away  from  0.  In  [9] , 
e is  then  set  equal  to  0 in  A - tB,  and  the  determinant  of  that  matrix  is  then  shown 
to  have  m sign  changes  as  r decreases  from  u.  Thus  the  conjectured  value  v bounds  all 
eigenvalues  of  A - XB  for  all  values  of  x.  The  upper  bounds  thus  obtained  are  v = .69 
for  to  = 3 and  v — .9  for  m — 4.  We  emphasize  that  the  B-splines,  matrix  entries,  and 
determinants  all  are  calculated  exactly,  using  the  Maple  ([6,  7])  computer  algebra  system, 
so  the  bounding  property  of  v is  rigorously  proven.  Since  the  bounds  we  obtain  apply  to 
the  spaces  Gi  and  B,  associated  with  any  one  of  the  a q,  they  satisfy  the  hypotheses  of 
Theorem  2.4  which  now  provides  our  conclusions.  □ 

Our  main  result  now  follows. 

Theorem  3.2  Let  the  hypotheses  be  those  of  Lemma  3.1.  In  addition,  let  P be  the 
orthogonal  projector  onto  the  space  of  n-th  order  real-valued  T -periodic  trigonometric 
polynomials,  where  n > 3N.  If  m = 3,  we  have  ||P+||2  < 2.4,  while  if  m = 4,  we  have 
||P+||2<4.5. 

Proof:  The  space  Q in  Lemma  3.1  consists  of  periodic  splines  with  uniformly  spaced 
knots.  Theorem  3.1  of  [10]  implies  that  ||/  - P||2  < (a/(  1 + a))1/2  where 

OO 

a=4]T(l/(l+2r))2''’. 

T — l 

In  [9]  we  use  this  formula  to  get  upper  bounds  of  .076  when  m = 3 and  .025  when  m = 4 
for  ||J  - P||2.  Taking  these  bounds  as  r)0  in  Theorem  2.3  and  taking  the  bounds  from 
Lemma  3.1  as  r)i  in  Theorem  2.3,  we  obtain  from  that  theorem  bounds  for  ||J  — P||2  of 
.907  for  m = 3 and  .974  for  m = 4.  Theorem  2.2  now  applies  to  produce  the  present 
results.  O 

Above,  we  required  n > 3 N-,  under  this  condition  we  can  get  our  simplest  and  most 
comprehensive  results.  Since  we  contemplate  applying  our  results  where  the  number  n 
of  useful  coefficients  may  be  limited,  we  have  tried  to  get  versions  of  Theorem  3.2  where 
n is  smaller  compared  with  N.  We  have  no  useful  versions  for  n < 3 N and  m = 4 (cubic 
splines).  The  following  result  for  quadratic  splines  may  be  useful.  To  formulate  it,  let 
t\  = max{j.T  - kx\N/T}.  In  the  previous  results,  the  separation  of  the  values  x from 
their  nearest  uniform  mesh  points  kx  was  unrestricted,  which  corresponds  to  = 1/2. 
Here,  we  can  get  results  for  quadratic  splines,  and  n > 2N,  provided  the  x are  more 
restricted;  our  methods  of  analysis  “blow  up”  for  n > 2 N as  fj  approaches  a number 
slightly  larger  than  .25. 
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Theorem  3.3  Let  m = 3 (quadratic  splines);  let  n > 2N.  Otherwise,  let  the  hypotheses 
he  those  of  Theorem  3.2.  Corresponding  to  the  list  0,  .1,  .2,  .25  for  values  of  e we  have 
the  list  of  values  1.7,2.1,3.9,16  as  bounds  for  j|P+||. 

Proof:  For  each  of  the  cases  for  e,  an  argument  similar  to  the  proof  of  Lemma  3.1 
applies  to  produce  a bound  rj i for  |j/  — G||2  where  Q now  is  defined  using  the  uniform 
knot  spacing  1/(2 N)  rather  than  1/(377).  The  only  difference  in  the  argument  is  that 
here,  a discontinuity  location  x always  stays  in  the  interval  \kx,  kx  + e\  h)  where  h = T/N, 
so  the  matrix  entries  and  determinants  can  be  treated  as  functions  of  e in  [0,  ei].  Each 
bound  rji  now  is  used  just  as  in  the  proof  of  Theorem  3.2,  to  get  the  present  bounds  for 

\\P+ 1|2-  □ 

4 Uniform  norm  bounds 

Using  representers  of  point  evaluation,  as  in  [8],  we  can  get  uniform  norm  bounds  for  P+, 
and  thus  for  A.  The  arguments  are  similar  to  those  in  [8].  The  main  difference  is  that 
there  the  mesh  is  uniform  and  the  order  m is  3.  The  constructions  of  representers  extend 
fairly  easily  to  the  present  case:  here  the  norms  of  representers  are  functions  both  of  the 
evaluation  point  and  the  location  of  the  discontinuity  nearest  to  the  evaluation  point. 
One  can  show  that  for  each  point  f€[0,T),a  spline  rt  exists  in  a space  U containing  Q, 
such  that  (rt,  q)  = q(t)  for  each  q € Q,  and  such  that  | |rt 1 12  < k / Vh  where  k = 5,  m — 3 
and  k = 7 ,m  — 4;  h — T/N  as  before.  The  computations  for  the  construction  and  bound 
calculations  are  in  [9] . Noting  that  \ff  / \fli  — \fN,  we  have 


P/lloc  < max||n||2p/||2  < (fc/x/piP+iuVril/lloo  < WiVllP+PI/IU 


When  N < 100  and  the  hypotheses  are  those  of  Lemma  3.1,  this  gives  P/||oo  < 
12011/Hoo  for  m = 3,  and  H-A/Hoo  < 315||/!|oo  for  m = 4. 

5 Example 


Fig.  1. 


f - Af,  no  noise 


f - Af,  1%  noise 


exact  f 


R.  K.  Wright 


We  illustrate  the  method  using  an  example  where  the  function  / is  27r-periodic  and  on 
[0, 27r)  consists  of  the  function  e-*/6  with  a piecewise  quadratic  added,  so  as  to  produce 
discontinuities  at  0,  .5, 1.5, 2.5,  and  4.  / is  a modification  of  an  example  in  [2];  for  conveni- 
ence we  have  shifted  that  example  left  by  1 unit,  and  we  have  added  the  exponential  term 
because  our  method  can  represent  a piecewise  quadratic  exactly  in  the  absence  of  noise. 
Exact  (up  to  17-decimal  digit  floating  point  error)  Fourier  coefficients  are  derived  from 
/ by  exact  integration  using  the  Maple  ([6,  7])  system.  Noisy  approximate  coefficients 
are  also  derived  by  sampling  / at  1024  equidistant  values  in  [0,  2tt],  adding  uniformly 
distributed  pseudo-random  noise  to  the  samples,  and  taking  the  discrete  Fourier  trans- 
form of  the  samples.  In  effect,  we  work  with  / + e where  e is  a perturbing  function. 
The  level  of  the  noise  is  set  so  that  the  discrete  I/2-norm  of  the  noise  vector  is  1%  of 
the  discrete  1,2-norm  of  the  vector  of  samples  of  /.  N = 45  and  thus  n = 135  are  the 
smallest  values  of  n and  N for  which  the  hypotheses  of  the  previous  section  are  satisfied. 
Using  these  values,  we  proceed  with  m = 4 (cubic  splines)  for  each  of  these  cases  for 
Fourier  coefficients.  Plots  of  / and  of  the  error  for  the  two  cases  appear  in  the  figure. 
The  ratio  ||/  - A(f  + e)| I2/I l/l I2  is  about  .005  for  the  case  of  1%  noise.  In  [9]  we  develop 
a probabilistic  estimate  of  .0037  for  the  ratio  of  | |e| I2/! |/| I2-  This  estimate  indicates  an 
Z/2-norm  noise  magnification  of  about  1.35-fold,  compared  with  the  upper  bound  of  4.5 
given  in  Theorem  3.2.  The  uniform  error,  for  noise-free  coefficients,  is  about  10“9;  com- 
putational experiments  show  this  is  dominated  by  truncation  error  in  approximating  the 
exponential  term.  In  [9]  we  do  the  corresponding  calculations  for  m = 3,  and  find  similar 
results  for  1%  noise,  with  larger,  but  still  small,  error  for  noise-free  coefficients. 

In  [9],  we  implement  Eckhoff’s  method  as  described  in  [3],  used  on  the  above  data. 
For  noiseless  data,  the  results  are  comparable  to  those  reported  by  Eckhoff  for  similar 
examples.  The  uniform  norm  error  seems  to  be  about  .06,  with  errors  at  jumps  somewhat 
smaller.  For  1%  noise,  the  results  of  Eckhoff’s  method  are  about  750-fold  in  error. 
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